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Abstract 

We establish a general vector field model to describe the role of transverse momentum currents 
in optical Magnus effect in the free space. As an analogy of mechanical Magnus effect, the circu- 
larly polarized wavepacket in our model acts as the rotating ball and its rotation direction depends 
on the polarization state. Based on the model we demonstrate the existence of a novel optical 
polarization-dependent Magnus effect which is significantly different from the conventional optical 
Magnus effect in that light-matter interaction is not required. Further, we reveal the relation be- 
tween transverse momentum currents and optical Magnus effect, and find that such a polarization- 
dependent rotation is unavoidable when the wavepacket possesses transverse momentum currents. 
The physics underlying this intriguing effect is the combined contributions of transverse spin and 
orbital currents. We predict that this novel effect may be observed experimentally even in the 
propagation direction. These findings provide further evidence for the optical Magnus effect in the 
free space. 
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I. INTRODUCTION 



The optical Magnus effect is the photonic version of the Magnus effect in classical me- 
chanical systems. As a result of the usual mechanical Magnus effect, a rotating ball falling 
through the air is deflected from the vertical to the direction of rotation. The influence of 
spin on the trajectory can be considered as the optical analogue of the mechanical Magnus 
effect. The spin photon to some extent can be regarded as a rotating ball, and then after 
multiple reflections during its propagation through a fiber the photon will be deflected from 
its initial trajectory jl, 2]. Hence this effect has received the name of optical Magnus effect 
or optical ping-pong effect. Recently, the physical nature of the optical Magnus effect is 
connected with spin-orbital interaction in a wave field, which essentially depends on the 
gradient of a refractive index of optical medium 

For beams propagating in the free space, it is generally believed that the optical magnus 
effect takes no place since the gradient of a refractive index is zero. In fact, the changes 
in transverse momentum currents can lead to the rotation phenomenon in paraxial light 
beams. The possible occurrence of a polarization-dependent rotation of beam centroid in 
the free space was already predicted by Borodavka and coworkers [4]. However, they do not 
connect such a rotation with the occurrence of nonzero transverse momentum currents. In 
our opinion, the optical Magnus effect manifests itself within vector field structure in the 
frame of classical electrodynamics. However, the Jones vector is not sufficient to describe 
the vectorial property of a finite beam, due to the longitudinal component [5[. Thus, it is 
necessary for us to establish a general vector field model to reveal the role of transverse 
momentum currents in optical Magnus effect. 

Meanwhile the same interaction also leads to other effects such as the spin Hall effect 

of light (SHEL) a a. 

The interesting effect has been recently observed in beam refrac- 
tion [8] and in scattering from dielectric spheres (9]. More recently, the SHEL was found to 
occur when a light beam is observed on the direction making an angle with the propagation 
axis [10| . This effect has a purely geometric nature and amounts to a polarization-dependent 
shift or split of the beam intensity distribution. As the tilting angle tends to zero, the split- 
ting effect vanishes. The possible reason for this is that the transverse intensity distribution 
of fundamental Gaussian beam is axially symmetric and the polarization-dependent rotation, 
if it exists, cannot be observed directly. We believe that such a polarization-dependent shift 
should be unavoidable when the beam possesses the transverse momentum currents. Thus, 
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our another motivation is to prove the intriguing optical Magnus effect can be observed even 
in the beam propagation direction. 

In this work, we want to explore what role of the transverse momentum currents play 
in the optical Magnus effect in the free space. First, we introduce the Whittaker scale 
potentials to describe the vector field structure for different polarization models. In the frame 
of classical electrodynamics, it is the rotating wavepacket but not the spin photon acting 
as the spin ball, which is significantly different from the previous works. Then, we uncover 
how the centroid of wavepacket evolves, and how the polarization state affects its rotation 
of the centroid trajectory. Finally, we examine what roles the spin and orbital currents 
play in the optical Magnus effect. A relation between transverse momentum currents and 
optical Magnus effect is obtained. We believe that our findings may provide insights into 
the fundamental properties of optical currents in beam propagation. 



II. VECTOR FIELD MODEL 



In order to reveal the role of the transverse momentum currents in optical Magnus effect, 
we begin to establish a general beam propagation model to describe the vector field structure. 
Figure [T] illustrates the geometry of the vector field structure in the Cartesian coordinate 
system. Following the standard procedure, the electric field E(r, u) is obtained by solving 
the vector Helmholtz equation 

(V 2 + A; 2 )E(r,u;) = 0, (1) 

where r = xe x + ye y + ze z and k = u/c is the wave number in the free space. The vector 
Helmholtz equation can be solved by employing the angular spectrum representation of the 
electric field as 

E(r) = J dk x dkyE(k x , k y ) 

x exp[i(k x x + k y y + k z z)]. (2) 

The transverse condition V ■ E(k x , k y ) = implies that the components of the angular spec- 
trum satisfy the relation k ■ E(fc x , k y ) = 0. In principle, due to the longitudinal component, 
the Jones vector is not sufficient to describe the vectorial properties of a finite beam. Hence 
it is necessary for us to define two mutually orthogonal vectors 




FIG. 1: (color online) Schematic illustrating the geometry of the vector field structure. The 
fixed unit vector u lies in the plane zox and makes an angle 6 with the propagation axis z. The 
propagation axis z perpendicular to the plane of xoy, and the wavevector k for an arbitrary angular 
spectrum component is perpendicular to the plane of I10I2 

describing the vector field structure. Here, the fixed unit vector u lies in the plane zox and 
makes an angle 9 with the propagation axis, and 



u = sin 9e x + cos 9e z 



(4) 



where — tt/2 < 9 < +tt/2. Note that this fixed unit vector is not a purely mathematical 



concept 
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121 ] . In fact, a perpendicular u to the pro pag ation axis corresponds to the 



uniformly polarized beam in the paraxial approximation [l3|, and a parallel u corresponds 
to the cylindrical vector beam [14j. The axis u that is neither perpendicular nor parallel to 
the propagation axis was observed by Hosten and Kwiat in a recent experiment |8|. 

In order to accurately describe the optical Magnus effect, we introduce the scaler Whit- 
taker potentials [5| to represent the vectorial field. Consequently, the angular spectrum can 
be decomposed along the two vectors and we have 



E(k x , k y ) = Vi(k x , k y )Ii + V 2 {k x , k y )h. 



(5) 
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Here Vj = (a (3) T f(k x ,k y ) (j = 1,2) are the scaler Whittaker potentials. The amplitude 
of the angular spectrum is referred to as 



f{k x ,k y ) = -7= exp 

V Z7T 



kl + ky 2 



(6) 



where Wq is beam-waist size. The matrix (a /3) T denotes the Jones vector which satisfies 
the normalization condition a/3* + = 1. The coefficients a and satisfy the relation 
a = i(a(3* — a* (3). The polarization operator o = ±1 corresponds to left and right circularly 
polarized light, respectively. It is well known that a circularly polarized beam can carry spin 
angular momentum ah per photon due to its polarization state 15 ]. 

From the viewpoint of Fourier optics [16], the Whittaker potentials Vj are given by the 
relation 

^'( r ) = ^2 J dk x dkyVj(k x ,ky) 

x exp[i(k x x + k y y + k z z)}. (7) 

It can be verified that both V\ and V2 satisfy the scalar Helmholtz equation: 

(V 2 + A; 2 )^(r) = 0. (8) 

On substituting Eq. (jSJ) into Eq. (j2J), the electric field can be expressed in terms of the 
Whittaker potentials 

E(r) = V x V x (uVx) - ikV x (uV 2 ). (9) 

In fact, after the angular spectrum on the plane z = is known, Eq. (Q together with 
Eqs. (Hj) and (jSJ) provides the expression of the electric field vectors as 

^ =sinS ^ + ^)v, 1 + C0 3^_ ic08 ^, (10 ) 

£ iwi = ^ e ^. + coa e^--i«nek^+ic 0a ek^, (11) 

y oxoy oyoz oz ox 

£ i wi — e ^ + ™ e {* + &h +i «* ek w- (i2) 

Here the superscript [W] represents the vectorial model given by the Whittaker potentials. 
We find that the electric components can be written as the Whittaker potentials and their 
first- and second-order derivatives. Up to now, we have established a general propagation 
model to describe the vector field structure. 
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It should be noted that the choice of the propagation models in the SHEL has highly 
debated a a. It is necessary for us to introduce the two polarization models, since op- 
tical Magnus effect shares the same physical mechanism with SHEL. Under the paraxial 
approximation of model [W], we can get the two different polarization models: 



E m 



oc 



a + /3y cot 9 J e x + ( f3 - cot 6 f e, 



ak x + (3k y 
k 



f{k x , k y ), 



(13) 



El- 



oc 



ae x + (3e y - 



ak x + f3k y 
k 



f(k x , k y 



(14) 



Evidently, model [II] can be obtained from model [I] under the condition 9 = it/ 2. In the 
context of wave optics, both model [I] and model [II] are exact solutions to the paraxial 
wave equation. However, the customary paraxial approximation has been shown to be 
incompatible with the exact Maxwell's equations [171 ] . The practical light beams consist of 
electromagnetic fields and hence are governed by Maxwell's equations. In recognition of this 
fact, it is sometimes suggested that the model [W] representing the Gaussian beam seems to 
be appropriate. From the experimental viewpoints, even an ideal polarizer will transmit part 
of a cross-polarized wave [l^Q], and thus cannot produce a pure linear polarization state. 
From the theoretical viewpoints, the linear polarization is a paraxial approximation solution 
of the Maxwell's equations. When we go beyond the paraxial approximation and consider 
the lowest order corrections, the field is elliptical polarization in the cross section 20, 121 ]. 



Thus, it is desirable to consider the mode [W] in the optical Magnus effect. As shown in the 
following, we will compare the results given by the three polarization models. 

The changes in the transverse momentum currents can be used to explain the physics 
behind the rotation phenomenon of the wavepacket, which results in a intensity redistribu- 
tion over the wavepacket cross section 22|, ]23j. The time-averaged linear momentum density 



associated with the electromagnetic field can be shown 



24) to be 



p[M]( r ) = J_ Re [ E [M]( r ) x H^r)], 



(15) 



Here M = W, I, II denotes different polarization models. The magnetic field can be obtained 
by H' M J = — ik~ l V x E' M l The momentum currents can be regarded as the combined 
contributions of spin and orbital parts 



P 



[M] 



(16) 



Here, the orbital term is determined by the macroscopic energy current with respect to an 
arbitrary reference point and does not depend on the polarization. The spin term, on the 
other hand, relates to the phase between orthogonal field components and is completely 
determined by the state of polarization {25]. In a monochromatic optical beam, the spin 
and orbital currents can be respectively written in the form 

pS ? M1 = Im[(E[ M ]-V)E[ M l*], (17) 

p£ I] =Im[E[ M ]*-(V)EM], (18) 

where E^* • (V)El M l = E [ ™ ] *VE [ ™ ] * + E [ ™ ] *VE [ ™ ] + E^WE^ is the invariant Berry 
notation [26] . It has been shown that both spin and orbital currents originate from the beam 
transverse inhomogeneity and their components are directly related to the azimuthal and 
radial derivatives of the beam profile parameters. However, the orbital currents are mainly 
produced by the phase gradient, while the spin currents are orthogonal to the intensity 



gradient [22j. As shown in the following, the spin and orbital currents will play different 
roles in the optical Magnus effect. 

The monochromatic beam can be formulated as a localized wavepacket whose spectrum is 
arbitrarily narrow j^. As a mechanical analogy, the circularly polarized wavepacket acts as 
the rotating ball in our model. To generate an asymmetric intensity distribution, we choose 
the angle of the fixed unit vector u as 6 = 7r/360. Note that such a vectorial beam should 



can be realized experimentally without technical difficulties [28|, |29j. In general, the rotation 
properties of wavepacket are expressed by the transverse momentum currents as shown in 
Fig. |2j Very surprisingly, the orbital currents are polarization-dependent in this polarization 
model [Figs. [2( a) and 12(d)] . This is due to the presence of polarization-dependent screw 
wavefront. For the left circular polarization a = +1, the total transverse momentum currents 
in the exterior part of wavepacket present an anticlockwise circulation, while in the inner 
part exhibit a clockwise circulation [Fig. 12(c)]. For the right circular polarization a = —1, 
the total transverse momentum currents present an opposite characteristics [Fig. [2(f)] . The 
inherent physics underlying this intriguing effect is the combined contributions of transverse 
spin and orbital currents. The wavepacket cannot be regarded as a rigid ball due to the 
opposite circulations. This is significantly different from the mechanical Magnus effect. 

For the model [W] under the condition 6 = n/2, which has also been verified by Hertz vec- 
tor method [30], the polarization-dependent rotation of the transverse momentum currents 
are plotted in Fig. [3] Very interestingly, the orbital currents are polarization-independent 




FIG. 2: (color online) The momentum currents for the model [W] under the condition 6 = 7r/360. 
First row: Left circularly polarized wavepacket a = +1. Second row: Right circularly polarized 
wavepacket a = — 1. The background distribution is depicted as the longitudinal currents and the 
green arrows is described as the transverse currents, [(a), (d)] Orbital momentum currents pj^. 
[(b), (e)] Spin momentum currents p^- [(c), (f)] Total momentum currents Pq^ + p<^. The 
cross section is chosen as z = Zr and the intensity is plotted in normalized units. 

in this polarization model [Figs. E(b) and [3(e)] . However, the transverse spin currents 
are polarization-dependent and the longitudinal spin currents are absent [Figs. [3(b) and 
[3(e)]. For the left circular polarization a = +1, the total transverse momentum currents 
present an anticlockwise circulation [Fig. [3(c)] . For the right circular polarization a = — 1, 
the total transverse momentum currents present a clockwise circulation [Fig. [3(f)]. The 
transverse intensity distribution of this model is axialiy symmetric and the polarization- 
dependent rotation, if it exists, cannot be observed directly. This is a possible reason why 
the polarization-dependent split is observed only on the plane which is not perpendicular to 
the propagation direction of the beam 10]. 
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x/A, x/A, x/A. 




x/A, x/A, x/A, 

FIG. 3: (color online) The momentum currents for the model [W] under the condition 9 = tt/2. 
First row: Left circularly polarized wavepacket a = +1. Second row: Right circularly polarized 
wavepacket a = — 1. The background distribution is depicted as the longitudinal currents and the 
green arrows is described as the transverse currents, [(a), (d)] Orbital momentum currents Pq^. 
[(b), (e)] Spin momentum currents p^- [(c), (f)] Total momentum currents +p^. The 
cross section is chosen as z = zr and the intensity is plotted in normalized units. 

III. OPTICAL MAGNUS EFFECT 



To illustrate the polarization-dependent rotation effect, we now determine the shift of 
wavepacket centroid, which is given by (r)t M l = (x)^e x + (y)^e y , with 



r)[M] = '77- . ( i 9 ) 



/ / rpf 1 \x,y,z)dxdy 
J J Pz(x,y, z)dxdy 



First let us examine the polarization-dependent rotation in the model [W] . By substituting 
Eqs. (ITU]) . ([Til and (fi~5|) into Eq. (03), we obtain the shifts as 



. . rwl z sin 26 , . 

( X ' 1 = - f 2« Mh ■ 2~~ff\ ' ( 20 ) 

2(cos^ + kz R sin 0) 

{y) m = ""7 M 3 ., , (2i) 

2(cos J + kz R sin 0) 

where Zr = kw\]/2 is the Rayleigh length. We find that the wavepacket centroid shifts 
a distance away from the propagation axis in xoy plane. The shift (x)' w ' is polarization 
independent, while the shift (y)^ depends on the polarization state a. Note that the shift 
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(x)[ w l can be regarded as a small angle inclining from the propagation axis 3lL |32| . However, 
the shift (|/)' w ' does not change on propagation. 

The angular positions of the wavepacket centroid indicates the rotation angle tan<^ w J = 



(x)^ I (y)^ w \ which is significantly different from the previous definition 33]. According to 
Eq. (fT9|) the rotation angle is given by 



J w] = arctan— . (22) 

OZr 

Very surprisingly, the rotation angle is proportional to the Gouy phase for circular polar- 
ization a — ±1. In all cross sections denoted by the distance from the z = plane, the 
transverse patterns are quite similar. However the overall scale changes and the pattern 
rotates as a whole. The former is caused by diffraction, and the latter is seen due to the 
existence of transverse momentum currents. 

Since the rotating wavepacket passes the distance between different cross-sections in time 
t = z/c, the instant angular velocity can be given by fl^(z) = dcp^/dt = cdcp^/dz, so 
that 

^(z) = c^— 2 . (23) 
a 2 z R + z l 

The rotation velocity of the wavepacket centroid decreases as the propagation distance in- 
creases. When the condition o = is satisfied and the rotation characteristics of wavepacket 
centroid vanish, i.e., Q^(z) = 0. This is why the linear polarized wavepacket cannot present 
the optical Magnus effect. However the linear polarization can be represented as a superpo- 
sition of two circularly polarized components (7|. As the result, the polarization-dependent 
split of the wavepacket intensity distribution arises. Thus, the same mechanism also leads 
to other effects such as the SHEL. 

Figure H] shows the polarization-dependent rotations of wavepacket centroid. At the 
plane z = const., the intensity of the wavepacket can be regarded as I^(r) = c 2 pi W ^(r). 
The wavepacket centroid (green spots) whose angular positions indicate the rotation angle 
if. For the left circular polarization o = +1, the wavepacket centroid exhibits a clockwise ro- 
tation [Figs. H[a)-|llc)]. For the right circular polarization a = —1, however, the wavepacket 
centroid presents an anticlockwise rotation [Figs. S^d)-Sl^f )] . During the wavepacket prop- 
agation from z = to z = +Zr, the rotation angle amounts to <£>' w l = an/ '4. The physics 
underlying this intriguing effect is the combined contributions of transverse spin and orbital 
currents. The novel polarization-dependent rotations differs from the conventional optical 
Magnus effect, in that light-matter interaction is not required. 
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FIG. 4: (color online) The polarization-dependent rotations of wavepacket centroid (green point) 
for the model [W]. First row: The wavepacket centroid undergoes a clockwise rotation for left 
circular polarization a = +1. Second row: The wavepacket centroid presents an anticlockwise 
rotation for right circular polarization a = — 1. Images are labeled by corresponding values of 
propagation length z and rotation angle . The intensity distributions are plotted in the plane 
[(a), (d)] z = 0, [(b), (e)] z = 0.5zr, and [(c), (f)] z = zr. Other parameters are the same as those 
in Fig. [21 

We are currently investigating the polarization mode [I]. On substituting Eqs. (j2D and 
(USD into Eq. flUD we have 

a cot 9 



<*> [I] =0 (y) 



(24) 



There exists an inherent transverse shift on propagation. This result coincides with that 
obtained by Li Uj with different methods. According to Eq. ( 1241) the rotation angle is given 
by 



[i] 



7T 

cr arctan — . 



(25) 



In this case, the rotation characteristics of wavepacket centroid vanish, i.e., Q^(z) = 0. 

We are now in a position to consider the polarization model [II] . On substituting Eqs. 
and (I14p into Eq. ( 1T91) we can determine 



(x)M=0, (y)M=0. 



(26) 



A further important point should be noted is that Eq. ([26]) can be obtained from Eq. (|2 
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FIG. 5: (color online) The polarization-dependent rotations of transverse momentum currents for 
the model [II]. First row: Left circularly polarized wavepacket a = +1. Second row: Right circu- 
larly polarized wavepacket a = — 1. The background distribution is depicted as the longitudinal 
currents and the green arrows is described as the transverse currents, [(a), (d)] Orbital momen- 
tum currents Pq • [(b), (e)] Spin momentum currents ply . [(c), (f)] Total momentum currents 
Pq 1 ' + . The cross section is chosen as z = zr and the intensity is plotted in normalized units. 



under the condition = n/2. In this model, the polarization-dependent rotation of the 
centroid cannot be observed directly. Thus, we attempt to explore an alternative way 
to describe the polarization-dependent rotation effect. When the momentum currents are 
included, the change rate of azimuthal angle with z axis is written as 34] 



dz 



pip 



rp z 



fill ' 



(27) 



Here, p v describes the momentum current that circulates around the propagation axis, and 
p z describes the momentum current that propagates along the +z axis. The result of Eq. f|T5|) 
can be expressed in term of azimuthal component, defined by 



Pf ] = -Px I] Sm f + Py l] COS <P, 



(28) 



By substituting Eqs. and (f5S]) into Eq. (|2"TI) and carrying out the integration, we obtain 



[ii] 



a arctan ■ 



zr 



(29) 
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This result is slightly different from the rotation of the centroid in model [W]. The instant 
angular velocity can be given by fi[ n l(z) = cd(p^/dz, so that 



= =ra- (30) 

R ' 

The rotation velocity of the momentum current decreases as the propagation distance in- 
creases. When the condition o = is satisfied, the rotation characteristics vanish. 

To obtain a clear physical picture, the polarization-dependent rotations of the momentum 
currents are plotted in Fig. El We find that the orbital currents are polarization-independent 
in this polarization model [Figs. EI a) and Eld)]. However, the spin currents are polarization- 
dependent [Figs.E^b) andE^e)]. For the left circular polarization a = +1, the total transverse 
momentum currents present an anticlockwise circulation [Fig. [51(c)] . For the right circular 
polarization a — — 1, the total transverse momentum currents present a clockwise circulation 
[Fig. E^f)]- The polarization-dependent rotations of momentum currents may provide an 
alternative way to illustrate the optical Magnus effect, although whether the momentum 
currents in the free space propagating along a curvilinear trajectory is still debated 
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Comparing Fig. [3] with Fig. El shows that there are no notable difference between the model 
[W] (under the condition 9 = 7r/2) and the model [II]. In fact, the latter can be regarded as 
the paraxial approximation of the former. 

IV. SPIN AND ORBITAL ANGULAR MOMENTA 

In the frame of classical electrodynamics, it is the circularly polarized wavepacket but 
not the spin photon acting as the rotating ball. Now a question arises: What roles the spin 
and orbital angular momenta play in the optical Magnus effect? We now in the position 
to analysis the angular momentum density for each of individual wavepacket, which can be 



written as |24j 



j [M| (r) = rxpM(r). (31) 
Within the paraxial approximation, the angular momentum can be divided into the spin 

.[M 

To 



and orbital angular parts j' M ^ = y^ 1 



jg 37], it follows that 




jg" = r X pS", 


(32) 


jf = r x pf. 


(33) 



This separation should hold beyond the paraxial approximation 
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We first consider the longitudinal angular momentum density j z which can be regarded 
as the combined contributions of spin and orbital parts: 

jte=Wcfi-yPc&> ( 34 ) 

Jsl ] = x p [ sS -ypfJ- (35) 

The longitudinal angular momentum will provide a simple way to understand why the 
circularly polarized wavepacket exhibits the optical Magnus effect. 

Figure E] presents the distribution of longitudinal angular momentum density for the 
model [W] . Very surprisingly, this polarization model of wavepacket possesses polarization- 
dependent orbital angular momentum density [Figs. |6^a) and Eld)]. However, the spin 
momentum density exhibits a significantly different distribution [Figs. Mjo) andE|e)]. For 
the left circular polarization a = +1, the total angular momentum density in the exterior 
part of wavepacket is positive j z > 0, while in the inner part is negative j z < [Fig. [6(c)] . 
This means that the outer part of the packet presents an anticlockwise rotation, the inner 
part undergoes a clockwise one. For the right circular polarization a = — 1, the total angular 
momentum density of wavepacket presents an oppositive distribution [Fig. Eff)]. Thus, 
the rotating wavepacket cannot be regarded as a rigid ball as in the mechanical Magnus 
effect. This is the reason why we choose the centroid as the reference point to describe the 
polarization-dependent trajectories. 

For the comparison, we plot the longitudinal momentum density of the model [II] in 
Fig. [3 By comparing with the model [W], we find that the orbital angular momentum 
density vanishes in the present polarization model [Figs. [3(a) and [3(d)]. Thus, only the spin 
momentum density exhibits a polarization-dependent distribution [Figs. [3(b) and [3(e)]. For 
the left circular polarization a — +1, the total angular momentum density first increases 
then decreases with the increase of r as shown in Fig. [3(c). For the right circular polarization 
a = —1, the total angular momentum density presents an oppositive distribution as shown 
in Fig. [7Jf). The transverse intensity distribution is axially symmetric and the polarization- 
dependent rotation, if it exists, cannot be observed directly. 

Now we want to enquire what roles the angular momentum play in the optical Magnus 
effect. To answer this question needs to discuss the transverse angular momentum. The time- 
averaged linear momentum and angular momentum, which can be obtained by integrating 
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FIG. 6: (color online) The distribution of the longitudinal angular momentum density for the 
model [W]. First row: Left circular polarization a = +1. Second row: Right circular polarization 
a = — 1. [(a), (d)] Orbital angular momentum density Jq^} ■ [(b), (e)] Spin angular momentum 
density Jg^ ■ [(c), (f)] Total angular momentum density Jq 1 } +J<^- The cross section is chosen as 
z = and the intensity is plotted in normalized units. Other parameters are the same as those in 

Fig. m 



over the whole x — y plane 24 ] 




p[M] = / / p [M] (r)dxd2/! (gg) 



j[M ] = y j j[M] (r)d;rd2/ _ (37) 

The transverse angular momentum components are given by 

jM = (y)pM-zPM, (38) 

J[M] = z p[M) _ {x)p [M]^ (39) 

In the z = plane, we have J^/pJ^ = (y)^ and J [ ™ ] /Pl M] = -(x)^. Thus, the 

transverse angular momentum can be obtained by measuring the position of wavepacket 

centroid [l0|. In order to reveal the rotation characteristics, it is necessary for us to know 

the transverse angular momentum in any cross section. We first consider the model [W], 

the transverse angular momenta are given by 

m _ no sin 20 [W] _ nzsm29 

Jx ~ 2kH R ' y ~ 2kH\ ' [ ' 
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FIG. 7: (color online) The distribution of the longitudinal angular momentum density for the model 
[II]. First row: Left circular polarization a = +1. Second row: Right circular polarization a = —1. 
[(a) and (d)] Orbital angular momentum density Jq^. [(b) and (e)] Spin angular momentum density 
3Sz' [( c ) an< ^ CO! Total angular momentum density jj^ • The cross section is chosen as z = 

and the intensity is plotted in normalized units. Other parameters are the same as those in Fig. 

In this case, the polarization-dependent rotation of the centroid is unavoidable since the 
wavepacket possesses transverse angular momentum. For the model [I] , we obtain 

ira cot 6 



J r . 



J [I] 



0. 



(41) 



k 2 z R 

This is the reason why the model [I] only exhibits a transverse shift. For the model [II], 
however 

4"] = o, Jf = 0. (42) 

Thus, the wavepacket centroid no longer presents the polarization-dependent rotation. The 
physics underlying this phenomenon is the absence of the transverse angular momentum. It 
should be noted that the SHEL can be noticeably enhanced when the wavepacket carries 



orbital angular momentum [39l44ll]. Further work is needed to uncover the optical Magnus 
effect of such a wavepacket in the free space. 

It should be mentioned that the recent advent of negative index metamaterials, also 



known as left-handed materials (LHMs) 



42 



can induce a reversed polarization- 



dependent rotation of the trajectory of the wavepacket centroid. Because of the negative 



index, we can expect a negative Rayleigh length in LHMs 

16 



44 



45]. It will be interesting for 



us to describe in detail how the wavepacket trajectory evolves in the LHMs. Recently, the 
technique of transformation optics has emerged as a means of designing metamaterials that 
can bring about unprecedented control of electromagnetic fields [46!]. It is possible that the 
trajectories of circularly polarized wavepacket can be controlled by introducing a prescribed 
spatial variation in the constitutive parameters. 

V. CONCLUSIONS 

In conclusion, we have established a general vector field model to describe the role of 
transverse momentum currents in optical Magnus effect in the free space. We have demon- 
strated the existence of a novel optical polarization-dependent Magnus effect which differs 
from conventional optical Magnus effect in that light-matter interaction is not required. In 
the optical Magnus effect, the circularly polarized wavepacket acts as the rotating ball, but 
is not identical to a rigid body. This is because different parts of the wavepacket present 
diverse rotation characteristics. For a certain circularly polarized wavepacket, whether the 
rotation is clockwise or anticlockwise depends on the polarization state. Such a polarization- 
dependent rotation is unavoidable when the wavepacket possesses transverse momentum 
currents. We predict that this novel effect may be observed experimentally even in the 
propagation direction. Our findings provide further evidence for the optical Magnus effect 
in the free space. Because of the close similarity in atom physics, condensed matter, and 
optical physics, we believe that the Magnus effect is not limited to electromagnetic fields, 
but extends to other research areas, such as atom, ion, and electron beams. 
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